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1.  INTRODUCTION 


The  reduced  cohesion  of  grain  boundaries  (GBs)  is  known  to  be  often  the  controlling  factor  limiting 
ductility,  and,  actually  the  performance  and  reliability  of  high-strength  metallic  alloys  (Briant  and  Baneiji 
1983;  Guttmann  and  McLean  1979).  Intergranular  embrittlement  in  metals  is  usually  caused  by  impurities 
segregating  towards  the  GBs  (Lee  et  al.  1984;  Meyers  et  al.  1965;  Seah  1980;  Seah  and  Hondros  1983). 
Impurities  present  in  bulk  concentrations  of  10"^-10“^  atomic  percent  can  result  in  a  dramatic  decrease 
of  plasticity,  drastically  degrading  mechanical  properties  of  metallic  alloys  and,  thus,  posing  significant 
technological  problems.  This  detrimental  effect  of  impurities  on  the  order  of  parts  per  million  may  be 
readily  understood:  a  simple  estimate  shows  that  a  ppm  amount  of  impurity  is  sufficient  for  saturating 
aU  the  GBs  in  a  typical  grain-size  polycrystal  (see  Table  1). 

Table  1.  Amounts  of  Impurity  Atoms  Sufficient  for  Forming  a  Monolayer  in  a  Grain  Boundary 


Grain  size  (mm) 

5 

20  50  100 

Amount  of  impurity  (atomic  ppm) 

60 

30 

15  6  3 

Recent  progress  in  developing  efficient  methods  of  first-principles  calculations  and  computational 
algorithms  made  possible  systematic  studies  of  the  role  of  impurities  in  intergranular  cohesion  on  the 
electron-ion  level.  Calculations  on  the  supercell  models  of  GBs  with  impurities  in  iron  (Wu,  Freeman, 
and  Olson  1992, 1993,  1994a,  1994b;  Krasko  and  Olson  1990, 1991;  Krasko  1992)  and  tungsten  (Krasko 
1993,  1994)  have  provided  an  in-depth  insight  into  mechanisms  of  GB  cohesion/decohesion  processes. 

Since  the  first-principles  electronic  calculations  on  low-symmetry  systems  (such  as  lattice  defects  or 
GBs)  are  still  extremely  complicated  and  costly,  semi-empirical  methods  based  on  solid  first-principles 
foundations  have  been  developed.  Among  them,  the  most  popular  is  the  embedded  atom  method  (EAM) 
(Daw  1989;  Daw  and  Baskes  1984;  Finnis  and  Sinclair  1984, 1986;  Daw  and  Baskes  1983,  1987).  This 
method  has  been  successfully  used  in  a  wide  variety  of  calculations. 

The  purpose  of  this  paper  is  to  elucidate  energetics  of  impurities  on  a  tantalum  GB  and  analyze  the 
effect  of  impurities  on  the  intergranular  cohesion  in  tantalum. 
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Rather  than  doing  sophisticated  first-principles  calculations  on  multiatom  models  of  a  GB,  we  have 
chosen  to  calculate  the  quantity  that  may  be  called  environment-sensitive  embedding  energy  (EE),  the 
energy  of  an  impurity  atom  in  an  atomic  environment  typical  for  a  GB  (capped-trigonal  prism). 
Knowledge  of  these  energies  for  different  impurities  enables  one  to  compare  the  relative  stability  of  a 
particular  impurity  in  the  tantalum  GB.  Our  approach  is  obviously  an  extension  of  the  "Effective 
Medium"  (EM)  theory  (Nprskov  and  Lang  1980;  Puska,  Nieminen,  and  Manninen  1981;  Jacobsen, 
Nprskov  and  Puska  1987;  Stott  and  Zaremba  1980, 1982;  Nprskov  1982).  Earlier,  the  EEs  of  impurities 
in  iron  (Krasko  and  Olson  1990,  1991;  Krasko  1992)  and  tungsten  (Krasko  1993, 1994)  were  calculated 
and  used  in  analysis  of  GB  stability  in  those  metals. 

Having  calculated  the  EEs  for  a  number  of  impurity  atoms,  one  can  use  this  information  on  the 
impurity  energetics  in  a  modified  EAM  approach  for  calculating  the  GB  relaxation.  The  latter  calculation 
enables  one  to  draw  important  conclusions  regarding  the  intergranular  cohesion  in  tantalum  in  the  presence 
of  a  definite  impurity  in  the  GB. 

2.  RESULTS  AND  DISCUSSION 

2. 1  Environment-Sensitive  EEs.  The  basic  idea  of  the  EM  approach  was  to  replace  the  low-symmetry 
system  consisting  of  an  atom  plus  a  host  matrix  by  a  high  symmetry  effective  system  of  the  atom  and  the 
homogeneous  electron  gas  of  a  density  equal  to  that  seen  by  the  atom.  The  energy  of  interest  was  called 
the  EE  and  was  equal  to  the  energy  difference  between  the  atom  embedded  in  the  electron  gas  and 
separately  the  isolated  atom  and  the  electron  gas. 

The  EM  theory,  as  a  first  and  very  crude  approximation,  completely  neglected  any  covalent  effects, 
though  one  could  expect  that  a  metalloid  impurity  in  a  transition  metal  would  develop  strong  sp-d 
hybridization,  resulting  in  covalent  bonds.  The  introduction  of  covalent  effects  via  perturbation  theory 
(Nprskov  1982),  resulted  in  sigitificant  corrections  in  the  embedding  fimction  for  hydrogen.  Further 
attempts  to  improve  upon  the  EM  method  were  undertaken  in  recent  years  (Raeker  and  DePristo  1990; 
Chetty  et  al.  1992). 

Improvements  in  the  EM  theory  actually  made  calculations  more  sophisticated,  spoiling  the  elegant 
simplicity  of  the  original  method.  Rather  than  introducing  corrections  to  the  EM  concept,  we  have  chosen 
to  perform  first-principles  calculations  on  a  simplified  model  of  a  GB  environment,  hydrostatically  varying 
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aU  the  characteristic  volumes,  thus  generating  a  series  of  "environment-sensitive"  EEs  as  a  function  of  the 
electron  charge  density  due  to  the  host  (tantalum)  atom  at  the  impurity  site.  Thus,  an  impunty  is  actually 
"embedded"  into  a  crystal  lattice  environment,  rather  than  into  an  electron  jeUium,  as  in  the  original  EM 
theory. 

The  model  chosen  for  the  GB  environment  is  an  eight-atom  hexagonal  supercell  (Ta^X,  where  X  is 
an  impurity  atom).  The  supercell,  together  with  the  capped-trigonal  prism  coordination  of  the  surrounding 
tantalum  atoms,  is  shown  in  Rgure  1. 


Figure  1.  The  Ta^X  hexagonal  sunenceU  emulating  a  typical  trigonal  prism  environment  of  tantalum  atoms 
in  the  ail)E3  GB:  a)  the  suoercell:  bl  ihe  trigonal-prism  coordination  (©  Ta.  •impurity). 

A  trigonal  prism  GB  configuration  foUows  from  the  theory  of  hard-sphere  packing.  Atomistic 
relaxation  studies  in  iron  have  also  shown  (Hashimoto  et  al.  1984a,  1984b)  that  an  impurity  atom,  such 
as  phosphoras  or  boron,  is  likely  to  occupy  an  interstitial  position  in  the  center  of  the  trigonal  prism 
formed  by  iron  atoms  in  the  GB  core  (even  if,  as  in  case  of  phosphorus  and  boron  in  iron,  the  impurity 
forms  a  substitutional  solid  solution  with  the  host).  The  hexagonal  superceh  has  been  chosen  both  because 
of  its  relatively  high  symmetry  and  its  emulation  of  a  (111)  13  GB  environment. 


3 


We  perfonned  the  spin-polarized,  scalar-relativistic,  Linear  Muffin  Tin  Oibitals  (LMTO)  (Skriver 
1984)  calculations;  the  Von  Barth  and  Hedin  (1972)  exchange-correlation  and  the  frozen-core 
approximations  were  also  used.  First,  a  series  of  calculations  (for  six  different  volumes)  was  performed 
with  an  impurity  absent  from  the  supercell  (i.e.,  an  empty  sphere  of  the  same  radius  as  that  of  the  radius 
of  the  impurity’s  Wigner-Seitz  sphere  was  substituted  for  the  latter).  Similar  calculations  were  then 
performed  for  each  of  the  impurities:  hydrogen,  boron,  carbon,  nitrogen,  oxygen,  phosphorus,  and  sulphur. 
The  EEs  were  defined  as  follows: 


EE  =  ECTag#)  -  ECTagO)  -  E(«),  (1) 

where  E(Tag#)  and  E(TagO)  are  respectively  the  energies  of  the  supercell  with  and  without  the  impurity 
(O  stands  for  an  empty  sphere  substituted  for  the  impurity  atom),  and  E(#)  is  the  energy  of  the  free 
impurity  atom.  In  order  to  make  the  calculations  more  consistent,  we  have  chosen  to  use,  as  E(#)s,  the 
values  of  E(Tag#)  -  E(TagO)  extrapolated  to  the  zero  charge  density  (n  =  0),  which  would  correspond  to 
the  energies  of  impurities  in  the  GB  environment  with  the  host  crystal  lattice  infinitely  expanded.  The 
EE  energies  (equation  [1])  as  a  function  of  n,  the  electron  charge  density  due  to  tantalum  atoms  at  the 
impurity  site,  are  presented  in  Figure  2. 
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Figure  2.  The  environment-sensitive  EEs  vs.  electron  charge  density  (in  atomic  imits  fa.u.1). 
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Plots  in  Figure  2  explain  the  phenomenon  experimentally  observed  in  many  metals,  the  "site- 
competition"  effect.  As  one  can  see,  in  the  range  of  electron  charge  density  typical  of  a  GB 
(0.015-0.025  a.u.),  carbon  has  the  lowest  energy  and,  thus,  would  push  the  other  impurities  off  the  GB. 
To  the  best  of  our  knowledge,  any  experimental  information  on  the  site-competition  effect  of  impurities 
in  tantalum  is  absent. 

The  plots  in  Figure  2  also  reveal  an  important  aspect  of  GB  impurity  behavior.  All  the  plots  (except 
for  that  for  hydrogen)  have  well-pronounced  minima.  The  positions  of  the  minima  correspond  to  electron 
density  at  the  impurity  site  due  to  the  surrounding  tantalum  atoms  that  would  occur  if  the  GB  were 
allowed  to  relax  in  such  a  way  as  to  minimize  the  impurity’s  energy.  The  minima  positions  systematically 
shift  toward  lower  densities  with  the  impurity  loosing  its  competitive  power.  Lower  charge  density  means 
a  more  "loose"  GB,  less  strong  and  more  prone  to  decohesion.  The  minimization  of  the  total  GB  energy 
(not  only  the  energy  of  the  impurity  atom)  gives  the  characteristic  charge  densities  which  are  somewhat 
higher  than  those  in  the  minima.  However,  from  this  point  of  view,  only  sulphur  and  phosphoras  are  the 
obvious  candidates  for  being  "decohesive,"  boron  and  oxygen  are  marginal,  while  carbon  and  nitrogen  may 
be  called  "cohesion  enhancers."  Calculations  of  GB  relaxation  shed  more  light  on  intergranular 
cohesion/decohesion  due  to  impurities. 

2.2  Grain  Boundary  Relaxation.  First-principles  calculations  are  still  too  time-consuming  and  costly 
to  be  used  for  the  investigation  of  GB  relaxation.  A  semiempirical  EAM  is  obviously  more  appropriate. 
In  order  to  find  both  the  embedding  functions  and  the  pair  potentials,  the  EAM  uses  experimental  data, 
such  as  cohesive  energies  and  elastic  moduli,  for  the  system  of  interest.  This  approach,  easily  applicable 
to  pure  metals  and  compounds,  may  become,  in  fact,  inaj^ropriate  if  the  effect  of  an  isolated  impurity 
atom  or  an  impurity  atom  on  a  GB  is  to  be  studied.  For  example,  one  can,  in  principle,  in  a 
metal-hydrogen  system,  use  the  embedding  function  and  the  pair  potentials  found  firom  experimental 
iirformation  on  that  metal’s  hydrides.  This  information,  however,  can  be  misleading  since  the  effea  of 
isolated  hydrogen  atoms  on  electronic  structure  and  cohesive  properties  of  the  metal  may  be  completely 
different  from  that  of  periodic  arrays  of  hydrogen  atoms  typical  of  hydrides.  In  Krasko  and  Olson  (1991), 
we  found  that  hydrogen  in  the  iron  GB  does  not  at  aU  contribute  its  electron  to  the  iron  d-valence  band, 
contrary  to  a  general  belief  that  in  transition  metals  the  hydrogen’s  electron  will  inevitably  go  to  a  d-band. 

In  order  to  resolve  this  difficulty,  we  have  chosen  to  calculate  the  energy  contributions  due  to  impurity 
atoms  in  the  GB  by  using  the  EEs  discussed  previously.  Since  the  EAM  functions  are  also  fundamentally 
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dependent  on  the  electron  charge  density  at  an  atom  site,  the  EEs  may  simply  be  added  to  the  EAM 
energy  of  the  host  atoms: 

E  =  Sr  V(R,R')  +  EE[n(Ri^p)].  (2) 

where  E^^j,  (n)  and  V(R,R’)  are  the  EAM  EE  and  the  pair  potential  as  fovmd  for  the  bulk  BCC  Ta  (we 
used  the  Finnis-Sinclair  functions  and  parameters  for  tantalum  [Finnis  and  Sinclair  1984,  1986]).  The 
third  term  is  die  energy  of  the  impurity  atom.  R  and  R'  are  the  positions  of  the  host  atoms,  Rjj„p  is  that 
of  the  impurity,  and  n(R)  and  nCR^j^p)  ate  the  electron  charge  densities  at  the  cite  of  a  host  atom  and  the 
impurity,  respectively.  The  electron  charge  density  at  a  given  site  can  be  taken  to  be  a  superposition  of 
the  finee  atom  charge  densities  or  found  from  more  sophisticated  procedures. 

As  mentioned  previously,  the  GB  environment  we  were  dealing  with  was  that  of  the  (1 1 1)  E3  tilt  GB. 
The  GB  structure  can  be  represented  as  a  succession  of  (111)  hexagonal  planes: 

....CBACBACBACBACBABCABCABCABCABC.... 

(the  GB  plane  is  marked  by  A).  The  CBABC  atomic  structure  of  the  core  of  the  GB — clean  (CL)  or  with 
an  impurity — ^was  just  emulated  by  the  superceU  and  is  shown  in  Figure  1.  hi  order  to  find  the  GB 
structure  corresponding  to  a  minimum  of  energy,  equation  (2),  the  interplanar  distances  were  varied,  while 
the  interatomic  spacings  and  the  structure  within  the  (111)  planes  were  unchanged. 

The  counterpart  of  the  nearest  neighbor  distance  in  the  bulk  BCC  lattice  (2.859  A)  is  the  distance 
between  two  tantalum  atoms  in  the  (111)  direction  Tal-Ta4  (see  Figure  1).  As  a  result  of  relaxation,  the 
CL  GB  this  distance  becomes  larger  than  in  die  bulk  (2.951  A),  while  the  shortest  distance  is  the  one 
between  two  Ta2  atoms  (Ta2-Ta2)  across  the  GB:  2.650  A.  In  the  CL  GB  there  is  a  significant  void 
(occupied  by  a  ★)  (see  Figure  3);  the  distance  between  Ta3  atoms  (Ta3-Ta3)  across  the  void  (and  the  GB 
plane)  is  quite  large:  3.504  A.  Thus,  in  the  CL  GB  the  strongest  interaction  is  Ta2-Ta2,  followed  by 
Tal-Tal  and  Tal-Ta2.  With  an  impurity  atom  in  place  of  ★,  the  interatomic  distances  and  interactions 
change  significantly.  Now  the  shortest  distance  is  that  between  an  impurity  atom  and  Ta3  (•-Ta3),  the 
distance  between  Ta2  atoms  across  the  GB  being  larger.  The  interaction  between  the  impurity  atom  and 
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Figures.  Schematic  arrangement  of  atomic  planes  in  GB  cores:  (a)  CL  GB;  (b)  GB  with  impurity. 

atom  TaS  becomes  of  utmost  importance.  It  is  actually  responsible  for  the  intergranular  cohesion;  this  is 
the  interaction  between  two  TaS  atoms  across  the  GB  via  the  impurity  atom  (see  Figure  S)  that  is  the 
strongest.  This  qualitative  picture  explains  why  even  the  tiniest  amounts  of  "bad"  impurities  at  the  GB 
may  cause  detrimental  worsening  in  metal’s  mechanical  properties.  On  the  other  hand,  an  impurity  with 
strong  interaction  with  the  host  atoms  will  result  in  cohesion  enhancement 

The  GB  relaxation  calculations  also  show  that,  as  in  the  case  of  the  iron  and  tungsten  GBs  (Krasko 
and  Olson  1990,  1991;  Krasko  1992,  1993,  1994),  the  interplanar  separations  oscillate  as  a  function  of 
distance  from  the  GB,  the  deformation  waves  decaying  by  the  10th-12th  plane  away  from  the  GB 
(Figure  3).  An  interesting  feature  of  the  CL  GB  relaxation  is  that  the  distance  between  the  second  and 
third  planes  is  less  than  half  of  the  (1 1 1)  interplanar  distance  in  bulk  BCC  tantalum  (0.430  A  vs.  0.954  A). 
This  "co-phase"  effect  (collapsing  of  two  (111)  planes)  has  been  also  found  in  iron  (Krasko  and  Olson 
1990,  1991;  Krasko  1992)  and  tungsten  (Krasko  1993,  1994).  The  simple  physical  explanation  of  it  is 
that  (as  one  can  see  from  Figures  1  and  3)  the  Ta3  atoms  across  the  GB  move  toward  each  other  simply 
to  decrease  the  volume  of  the  void  between  them.  When  an  impurity  atom  sits  in  the  void,  the  co-phase 
effect  disappears. 
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The  impurity  atoms,  cait)on  and  hydrogen,  cause  some  "damping"  of  the  relaxation  deformation  waves 
(i.e.,  decreasing  the  oscillation  amplitudes  with  respect  to  those  in  a  CL  GB).*  The  amplitude  increases 
with  oxygen,  phosphorus,  and  sulphur.  Damping  the  deformation  wave  may  be  interpreted  as  "cohesion 
enhancement,"  while  increasing  the  deformation  wave  oscillations  may  be  thought  of  as  resulting  in 
"decohesion." 

Carbon  and  nitrogen  are  known  to  be  GB  cohesion  enhancers  in  steels,  while  oxygen,  phosphoras,  and 
sulphur  (strong  embiitflers)  are  believed  to  weaken  the  GB  cohesioa  There  is  evidence  that  oxygen  in 
parts-per-million  concentrations  results  in  intergranular  embrittlement  of  tantalum  (Diaz  and  Reed-Hill 
1979).  At  the  same  time,  the  addition  of  carbon  seems  to  somewhat  neutralize  the  adverse  effect  of 
oxygen  (Kumar,  Mosheim,  and  Michaluk  1994).  The  latter  may  be  a  result  of  the  site-competition  eflfect 
(a  site-competition  cleansing).  Unfortunately,  reliable  experimental  information  on  effect  of  impurities  on 
tantalum  mechanical  properties  is  still  lacking. 

Hydrogen  is  one  of  the  worst  embrittlers  in  many  metals.  The  EAM  calculations  (Raeker  and  DePristo 
1990)  have  shown  that  a  hydrogen  atom  on  a  nickel  GB  does  weaken  the  metallic  bond  across  the  GB, 
lowering  the  fracture  stress  by  some  15%.  The  first-principles  calculations  in  iron  (Daw  and  Baskes  1983, 
1987)  suggest  that  one  of  the  decohesive  factors  is  a  weak  Fe-H-Fe  bond  across  the  GB.  Experimental 
information  on  the  effect  of  hydrogen  on  tantalum  mechanical  properties  is  also  lacking. 

2.3  Grain  Boundary  Stability.  From  a  thermodynamic  point  of  view  (Rice  and  Wang  1989; 
Anderson,  Wang,  and  Rice  1990),  the  irnpurity’s  embrittling  potency  depends  on  the  difference  between 
the  free  energies  of  the  impurity’s  segregation  on  the  initial  GB  and  on  the  two  free  surfaces  emerging 
upon  fracture.  The  higher  the  difference,  the  stronger  the  embrittling  potency  of  the  impurity.  As  a  less 
rigorous  but  simpler  criterion,  in  Seah  (1980)  and  Seah  and  Hondros  (1983)  the  sublimation  energy 
differences  between  the  host  and  impurity  were  compared  in  an  ideal  solution  model  for  over  60  elements. 
According  to  that  criterion,  in  tantalum,  only  nitrogen,  boron,  and  carbon  are  cohesion  enhancers  (with 
carbon  being  the  strongest)  (Figure  5). 


*  At  the  present  time,  we  cannot  explain  such  a  behavior  of  hydrogen;  its  EE  curve,  in  Figure  2,  is  quite  different  from  that  of 
carbon,  though  the  relaxed  interplanar  distances  are  almost  identical. 
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Figure  5.  The  diagram:  sublimation  energies  vs.  atom  size  calculated  from  an  ideal  solution  model  in 
Seah  (1980)  and  Seah  and  Hondros  (1983).  (Elements  above  and  below  the  broken  line, 
respectively,  increase  and  reduce  the  fracture  energy  in  tantalum.) 


The  effects  of  impurities  on  GB  stability  can  be  also  analyzed  by  simply  comparing  the  GB  energy 
differences  between  the  GB  with  impurities  and  the  clean  GB.  The  corresponding  values  for  the 
impurities  discussed  are  plotted  in  Figure  6.  One  can  see  that  the  AE  values  are  more  negative  again  for 
boron,  carbon,  and  nitrogen  than  can  be  thought  of  as  "cohesion  enhancers,"  and  less  negative  for  the 
potential  embrittlers;  hydrogen,  oxygen,  phosphorus,  and  sulphur.  Hydrogen  should  be  the  worst 
embrittler;  its  AE  is  only  slightly  negative  (a  positive  AE  would  mean  instability  of  the  GB  with  impurity 
at  absolute  zero  temperature). 
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Figure  6.  AE  =  EQ3(X)-EQg(cln),  the  energy  diiference  between  the  GB  with  impurity  X  and  clean  GB, 
vs.  die  periodic  chart  grout?  number. 

3.  CONCLUSIONS 

The  environment-sensitive  embedding  energies  (EEs) — ^the  energies  of  an  impurity  atom  in  an 
environment  typical  of  a  GB  in  tantalum  as  a  function  of  electron  charge  density  at  the  impurity 
site — ^were  calculated  for  a  number  of  impurities.  The  knowledge  of  these  energies  enabled  us  to  predict 
the  site  competition  effect  in  tantalum,  and  draw  conclusions  about  the  relative  stability  of  GBs  with 
different  impurities.  The  calculation  of  GB  relaxation  also  revealed  some  interesting  features  of  the  GB. 
Since  the  EEs  can  be  easily  calculated  for  various  host  environments,  a  modified  EAM  can  be  used  to 
analyze  the  impurity’s  energetics  in  a  variety  of  situations. 


11 


4.  REFERENCES 


Anderson,  P.  M.,  J-S.  Wang,  and  J.  R.  Rice.  "Thermodynamic  and  Mechanical  Models  of  Interfacial 
Embrittlement."  Irmovations  in  Ultrahigh-Streneth  Steel  Technology,  edited  by  G.  B.  Olson, 
M.  Azrin,  and  E.  S.  Wright,  34th  Sagamore  Army  Research  Conference  Proceedings,  p.  619, 1990. 

Briant,  C.  L.,  and  S.  K.  Baneiji  (eds.).  Embrittlement  of  Engineering  Alloys.  New  York:  Academic 
Press,  vol.  21,  1983. 

Chetty,  N.,  K.  Stokboro,  K.  W.  Jacobsen,  and  J.  K.  Nprskov.  "An  Ab  Initio  Potential  for  Solids." 
Physical  Review,  vol.  B46,  p.  3798,  1992. 

Daw,  M.  S.  "Model  of  Metallic  Cohesion:  The  Embedded  Atom  Method."  Rivsical  Review,  vol.  B39, 
p.  7441, 1989. 

Daw,  M.  S.,  and  M.  I.  Baskes.  "Semi-Empirical  Quantum-Mechanical  Calculations  of  Hydrogen 
Embrittlement  in  Metals."  Physical  Review  Letters,  vol.  50,  p.  1285,  1983. 

Daw,  M.  S.,  and  M.  I.  Baskes.  "Embedded  Atom  Method:  Derivation  and  Application  to  Impurities, 
Surfaces,  and  Other  Defects  in  Metals."  Physical  Review,  vol.  B29,  p.  6443,  1984. 

Daw,  M.  S.,  and  M.  I.  Baskes.  "Applications  of  the  Embedded  Atom  Method  to  Hydrogen 
Embrittlement."  Chemistry  and  Physics  of  Fracture,  edited  by  R.  H.  Jones  and  R.  M.  Latanision, 
Martinus  Nijhoff,  p.  196,  1987. 


Diaz,  A.  E.,  and  R.  E.  Reed-HUl.  "Evidence  for  Slow  Strain-Rate  Embrittlement  in  Tantalum  Due  to 
Oxygen."  Scripta  MetaUurgica.  vol.  13,  p.  491,  1979. 

Finnis,  M.  W.,  and  J.  E.  Sinclair.  "A  Simple  Empirical  N-Body  Potential  for  Transition  Metals." 
Philosophical  Magazine,  vol.  A50,  p.  45,  1984. 

Finnis,  M.  W.,  and  J.  E.  Sinclair.  "A  Simple  Empirical  N-Body  Potential  for  Transition  Metals." 
Philosophical  Magazine,  vol.  A53,  p.  161,  1986. 

Guttmarm,  M.,  and.  D.  McLean.  Interfacial  Segregations.  Edited  by  W.  C.  Johnson  and  J.  M.  Blakely, 
Metals  Park,  OH:  American  Society  for  Metals,  p.  261,  1979. 

Hashimoto,  M.,  Y.  Ishida,  R.  Yamamoto,  and  M.  Doyama.  "Atomistic  Studies  of  Grain  Boimdary 
Segregations  in  Fe-P  and  Fe-B  Alloys."  Acta  Metallurgica.  vol.  32,  p.  1,  1984a. 

Hashimoto,  M.,  Y.  Ishida,  R.  Yamamoto,  M.  Doyama,  and  T.  Fujiwara.  "Atomistic  and  Electronic 
Structures  of  a  Grain  Boundary  in  Iron  with  Impurity  Segregation."  Surface  Science,  vol.  144,  p.  182, 
1984b. 

Jacobsen,  K.  W.,  J.  K.  Nprskov,  and  M.  J.  Pudca.  "Interatomic  Interactions  in  the  Effective-Medium 
Theory."  Physical  Review,  vol.  B35,  p.  7423,  1987. 


12 


Krasko,  G.  L.  "Environment  Sensitive  Embedding  Energies  and  Grain  Boundary  Relaxation  in  Iron." 
Structuro  and  Properties  of  Interfaces  in  Materials,  edited  by  W.  A.  T.  Staik,  U.  Dahmen,  and 
C.  L.  Briant,  Mat  Res.  Soc.  Symp.  Proc.,  vol.  238,  p.  481,  Pittsburgh,  PA,  1992. 

Krasko,  G.  L.  "Site  Competition  Effect  of  Impurities  and  Grain  Boundary  Stability  in  Iron  and  Tungsten." 
Scripta  Metallurgica  et  Materialia,  vol.  28,  p.  1543,  1993. 

Krasko,  G.  L.  "Effect  of  Impurities  on  the  Electronic  Structure  of  Grain  Boxmdaries  and  Intergranular 
Cohesion  in  Timgsten."  International  Journal  Refractory  Metals  and  Hard  Materials,  vol.  12,  p.  251, 
1994. 

Krasko,  G.  L.,  and  G.  B.  Olson.  "Effect  of  Boron,  Carbon,  Phosphoras,  and  Sulphur  on  Intergranular 
Cohesion  in  Iron."  Solid  State  Communications,  vol.  76,  p.  247,  1990. 

Krasko,  G.  L.,  and  G.  B.  Olson.  "Effect  of  Hydrogen  on  the  Electronic  structure  of  a  Grain  Boundary  in 
Iron."  Solid  State  Commimications,  vol.  79,  p.  113,  1991. 

Kumar,  R.,  C.  E.  Mosheim,  and  C.  A.  Michaluk.  "Effect  of  IntermetalUc  Compounds  on  the  Properties 
of  Tantalum."  High  Temperature  Silicides  and  Refractory  Alloys,  edited  by  C.  L.  Briant, 
J.  J.  Petrovic,  B.  P.  Bewlay,  A.  K.  Vasudevan,  and  H.  A.  Lipsitt,  Materials  Research  Society  Symposia 
Proceedings,  vol.  322,  p.  413,  Pittsburgh,  PA,  1994. 

Lee,  D.  Y.,  E.  V.  Barrera,  J.  P.  Staik,  and  H.  L.  Marcus.  "The  Influence  of  Alloying  Elements  on 
Impurity  Induced  Grain  Boundary  EmbrittlemenL"  Metallurgical  Transactions  A,  vol.  15A,  p.  1415, 
1984. 

Meyers,  Jr.,  C.  L.,  G.  Y.  Onoda,  A.  V.  Levy,  and  R.  J.  Kotfila.  "Role  of  the  Grain  Boundaries  in  the 
Ductile-Brittle  Transition  Behavior  of  BCC  Refractory  Metals."  Transactions  of  Metallurgical  Society 
of  AIME.  vol.  233,  p.  720,  1965. 

Nprskov,  J.  K.  "Covalent  Effects  in  the  Effective-Medium  Theory  of  Chemical  Bonding:  Hydrogen  Heats 
of  Solution  in  the  3d  Metals."  Rivsical  Review,  vol.  B26,  p.  2875,  1982. 

Nprskov,  J.  K.,  and  N.  D.  Lang.  "Effective-Medium  Theory  of  Chemical  Bonding."  Physical  Review, 
vol.  B21,  p.  2136,  1980. 

Puska,  M.  J.,  R.  M.  Nieminen,  and  M.  Manninen.  "Atoms  Embedded  in  an  Electron  Gas:  Immersion 
Energies."  Physical  Review,  vol.  B24,  p.  3037, 1981. 

Raeker,  T.  J.,  and  A.  E.  DePristo.  "Corrected  Effective  Medium  Calculations  of  Chemisorption  of  H  and 
N  on  Fe  (100),  Fe  (110)  and  W  (110)."  Surface  Science,  vol.  235,  p.  84,  1990. 

Rice,  J.  R.,  and  J.-S.  Wang.  "Embrittlement  of  Interfaces  by  Solute  Segregations."  Materials  Science  and 
Engineering  A.  vol.  107,  p.  23, 1989. 

Seah,  M.  P.  "Adsorption-Induced  Interface  Decohesion."  Acta  Metallurgica,  vol,  8,  p.  955,  1980. 

Seah,  M.  P.,  and  E.  D.  Hondros.  "Atomistic  Mechanisms  of  Intergranular  EmbrittlemenL"  Atomistics 
of  Fracture,  edited  by  R.  M,  Latanision  and  J.  R.  Pickens,  New  York:  Plenum  Press,  p.  855,  1983. 


13 


Skriver,  H.  L.  The  LMTO  Method.  Berlin:  Springer,  1984. 


Stott,  M.  J.,  and  D.  M.  Zaremba.  "Quasiatoms:  An  Approach  to  Atoms  in  Nonuniform  Electronic 
Systems."  Rivsical  Review,  vol.  B22,  p.  1564,  1980. 

Stott,  M.  J.,  and  E.  Zaremba.  "Energies  of  Light  Atoms  in  Uniform  Electronic  Host  Systems."  Canadian 
Journal  of  Rivsics,  vol.  60,  p.  1145,  1982. 

Von  Barth,  U.,  and  L.  Hedin.  "A  Local  Exchange-Correlation  Potential  for  the  Spin-Polarized  Case." 
Journal  of  Rivsics.  vol.  C5,  p.  1629,  1972. 

Wu,  R.,  A.  J.  Freeman,  and  G.  B.  Olson.  "On  the  Electronic  Basis  of  the  PhosiAorus  Intergranular 
Embritflement  of  Iron."  Journal  of  Materials  Research,  vol.  7,  p.  2403,  1992. 

Wu,  R.,  A.  J.  Freeman,  and  G.  B.  Olson.  "Effect  of  P  and  B  Adsorption  on  Fe  (111)  Free  Surface." 
Physical  Review,  vol.  B47,  p.  6855,  1993. 

Wu,  R.,  A.  J.  Freeman,  and  G.  B.  Olson.  "Nature  of  Phosphorus  Intergranular  Embrittlement  for  the  Fe 
[110]  (111)  Grain  Boundary."  Physical  Review,  vol.  B50,  p.  75, 1994a. 

Wu,  R.,  A.  J.  Freeman,  and  G.  B.  Olson.  "First  Principles  Determination  of  Embrittlement  Behavior  of 
P  and  B  in  Fe  [110]  (111)  Grain  Boundary."  Science,  vol.  265,  p.  376,  1994b. 


14 


NO.  OF 

COPIES  ORGANIZATION 


2  DEFENSE  TECHNICAL  INFO  CTR 
ATTN  DTIC  DDA 

8725  JOHN  J  KINGMAN  RD 
STE0944 

FT  BELVOIR  VA  22060-6218 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTN  AMSRL  OP  SD  TA 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTN  AMSRL  OP  SD  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
ATTN  AMSRL  OP  SD  TP 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 


ABERDEFN  PROVING  GROUND 

5  DIR  USARL 

ATTN  AMSRL  OP  AP  L  (305) 


15 


NO.  OF 

COPIES  ORGANIZATION 

ABERDEEN  PROVING  GROUND 

30  DIR  USARL 

ATTN  AMSRL  MA  CC 
G  KRASKO 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  conunents/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number/Author  ARL-TR-1091  (Krasko) _ Date  of  Report  May  1996 - 

2.  Date  Report  Received  _ _ _ _ _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.)  _ _ _ _ _ _ — - - - 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) - 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ _ _ _ _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.)  - - - - - - - - - - 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


